Flux front penetration in disordered superconductors 



O 

o 

(N 



Stefano Zapperi 1 , Andre A. Moreira 2 and Jose S. Andrade Jr. 2 
1 INFM sezione di Roma 1, Dipartimento di Fisica, Universitd "La Sapienza" , P.le A. Moro 2, 00185 Roma, Italy. 
2 Departamento de Fisica, Universidade Federal do Ceard, 60451-970 Fortaleza, Ceard, Brazil. 

(February 1, 2008) 

We investigate flux front penetration in a disordered type II superconductor by molecular dy- 
namics (MD) simulations of interacting vortices and find scaling laws for the front position and 
the density profile. The scaling can be understood performing a coarse graining of the system and 
writing a disordered non-linear diffusion equation. Integrating numerically the equation, we observe 
a crossover from flat to fractal front penetration as the system parameters are varied. The value of 
the fractal dimension indicates that the invasion process is described by gradient percolation. 
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The magnetization properties of type II superconduc- 
tors have been studied for many years, but the interest 
in this problem has been renewed with the discovery of 
high temperature superconductors (!]]. The magnetiza- 
tion process is usually described in terms of the Bean 
model [[| and its generalizations: flux lines enter into 
the sample and, due to the presence of quenched disor- 
der, give rise to a steady flux gradient. While the Bean 
model provides a consistent picture of average magneti- 
zation properties, such as the hysteresis loop and thermal 
relaxation effects || , it does not account for local fluctu- 
ations in time and space. It has been recently observed 
that flux line dynamics is intermittent, taking place in 
avalanches and flux fronts are not smooth [||-|7|. In 
particular, it has been shown that the flux front crosses 
over from flat to fractal as a function of material param- 
eters and applied field (6). 

A widely used modeling strategy to describe the fluc- 
tuations around the Bean state consists in molecular dy- 
namics (MD) simulations of interacting flux lines, pinned 
by quenched random impurities IqTlq]. With this ap- 
proach it has been possible to reproduce flux profiles 
[ HI , hysteresis JlCj ] , avalanches ||,[l2|,[l3| and plastic flow 
[ s|jj3| . One of the aims of these studies jl0| is to estab- 
lish precise connections between the microscopic models 
and the macroscopic behavior, captured for instance by 
generalized Bean models. A different approach treats the 
problem at mesoscopic scale, describing the evolution of 
interacting coarse-grained units |i"4| , |l5| , supposed to rep- 
resent the system at an intermediate scale. While these 
models give a faithful representation of several features 
of the problem, the connection with the underlying mi- 
croscopic dynamics remains unexplored. 

In this letter, we investigate the invasion of magnetic 
flux into a disordered superconductor, a problem that 
has recently been the object of intense experimental re- 
search |5|-[7j . We first analyze the problem by MD simu- 
lations, in analogy with Ref. ifTof , and obtain scaling laws 
relating the front position and the flux profile to the pin- 
ning strength. In order to understand these results, we 
perform a coarse-graining of the equation of motion and 
obtain a non-linear diffusion equation for the flux line 



density. In the absence of pinning, the equation reduces 
to the model discussed in Refs. [jl(|[l7|. This model has 
been analytically [ fU3||i~7| ] solved to provide expressions for 
the initial dynamics of the front that are in agreement 
with our MD simulations. We show that when quenched 
disorder is included in the diffusion equation, the flux 
front roughens and is eventually pinned. Varying the pa- 
rameters of the model (applied field, disorder, interaction 
strength), the fluctuations of the front display a crossover 
from flat to fractal that is consistent with experimental 
observations |3|-[jJ. The value of the fractal dimension 
suggests that the strong disorder limit is described by 
percolation. 

In an infinitely long cylinder, flux lines can be mod- 
eled as a set of interacting particles performing an over- 
damped motion in a random pinning landscape [p|-p^[. 
The equation of motion for each flux line i can be writ- 
ten as 

Tin = An - fj) + E G[(Rp - n)/l] + ri(n,t), (l) 



where the effective viscosity can be expressed in terms 
of material parameters V = $0^2/ PnC 2 ■ Here, $0 is 
the magnetic quantum flux, c is the speed of light, p n 
is the resistivity of the normal phase and H c i is the up- 
per critical field. The first term on the right hand side 
represents the vortex-vortex interaction and it is given 
by J(r) = [$^/(87rA 3 )]Xi(|f|/A)f where the function K t 
is a Bessel function decaying exponentially for \r\ > A 
and A is the London penetration length jl^] . The second 
term on the right hand side accounts for the interaction 
between pinning centers, modeled as localized traps, and 
flux lines. Here, G is the force due to a pinning cen- 
ter located at R p , I is the range of the wells (typically 
I <C A), and p = l,...,N p (N„ is the total number of 
pinning centers). While Refs. P-pd| employed parabolic 
traps, we decided to avoid discontinuities in the force 
and used instead G(x) = — f x(\x\ — l) 2 , for \x\ < 1 and 
zero otherwise. Finally, we add to the model an uncorre- 
cted thermal noise term 77, with zero mean and variance 
(rf) = kbT/T. Although in the present simulations we 
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restrict ourselves to the case T = (see Ref. Jl8j for the 
implementation of thermal noise in MD simulations), we 
will discuss the effect of temperature in the coarse grained 
description of the dynamics. 

We perform MD simulations based on Eq. ([!]) and 
analyze the flux front propagation for different values 
of the pinning strength fo. We use N — 5000 flux 
lines in a system of size (L x = 800X,L y = 100A [pof), 
with N p = 800000 Poisson distributed pinning centers of 
width I = A/2, corresponding to a density of n — 10/A 2 . 
The injection of magnetic flux into the sample is imple- 
mented similarly to Ref. jl0| , concentrating at the begin- 
ning of the simulation all the flux lines in a small strip, 
parallel to the y direction, of length L' = 10~ 2 A and im- 
posing periodic boundary conditions in both directions. 
The front position x p is identified as the x coordinate 
of the most advanced particle in the system at differ- 
ent times. In Fig. 1 we show that x p grows initially 
with time as i 1 / 3 and eventually saturates to a value 
£ p which increases as the strength of the pinning cen- 
ters fo is decreased. The data collapse shown in the 
inset of Fig. 1 indicates that the front pinning length 

1/2 

£ p scales as f . When the front saturates, we mea- 
sure the density profiles E^] which can be collapsed us- 
ing the scaling form p(x) = f 1 ^ 2 F(xf^ 2 ) (see Fig. 2). 
In addition, we vary the density of pinning centers, us- 
ing UqX 2 — 2.5, 5, 7.5, 10, and find that similar scaling 
collapses hold if fo is replaced by fo^/no- 

In order to understand these results, we perform a 
coarse graining of Eq. ([!]) , starting from the Fokker-Plank 
equation for the probability distribution of the flux line 
coordinates P(r\, ....,fjv,t) 

BP 

r— =^V i (-/ i P + fc B TV i P), (2) 

i 

where fi is the force on the particle i given by Eq. (Q). 
Next, we introduce the single particle density p(r, t) = 
(y] { S 2 (r — fi)), where the average is done over the distri- 
bution P{f\ 1 ....,rjsr,t). The evolution of p can be directly 
obtained from Eq. (g) and is given by 

T % = ~ f ( / dVj V-Op (2) (r>"*',*) 

- G[(R P - r)//]p(f, t)) + k B TV 2 p, (3) 
p 

where p^(r,r',t) is the two-point density, whose evolu- 
tion depend on the three-point density and so on. The 
simplest truncation scheme involves the approximation 
p( 2 \r, r',t) ~ p(r, t)p(f ', t). We then coarse grain the 
equation considering length scales larger than A. This 
can be done expanding J in Fourier space, keeping only 
the lowest order term in q, and retransforming back in 
real space. The result reads 

J d 2 r'J(r~r')p(r',t)~~aVp(r,t), (4) 



where a = j d 2 rr ■ J(r)/2 = <I>o/4. 

The coarse graining of the disorder term is more sub- 
tle. A straightforward elimination of short wavelength 
modes would give rise, as in the previous case, to a 
random force F c (r) — — gVn, where n is the coarse 
grained version of the microscopic density of pinning cen- 
ters h(r) = J2 P 5 2 {r — R p ) and g cx f . This method can 
not be applied for short-range attractive pinning forces 
as the one we are investigating. In this case, short wave- 
length modes yield a macroscopic contribution to pinning 
that can not be neglected. Consider for instance the flow 
between two coarse grained regions: short-range micro- 
scopic pinning forces give rise to a macroscopic force that 
should always oppose the motion, while the random force 
derived above could in principle point in the direction of 
the flow. In other words, F c (r) should be considered as 
a friction force whose direction is always opposed to 
the driving force Fd (in our case Fd — aVp) and whose 
absolute value is given by |(?Vro| for \Fd\ > |<?Vn| and to 
\F d \ otherwise p|. 

Collecting all the terms, we finally obtain a disordered 
non-linear diffusion equation for the density of flux lines 

T^ = V(apVp-pF c ) + k B TV 2 p. (5) 
ot 

The boundary conditions representing our MD simula- 
tions correspond to a constant number of flux lines ML y 
injected into the system, so that the total density is con- 
served J dxdyp(x,y,t) = ML y . With these boundary 
conditions, the T = behavior depends on a single effec- 
tive coupling constant go = g^/no/ (aM) as can be shown 
rescaling Eq. (|5|) and the conservation law as 

p = p/M t = t/(TaM 2 ) n = ny^O- (6) 

In the limit go — and T = 0. Eq. (||) has been solved 
by Bryskin and Dorogotsev |l6| using similarity meth- 
ods and the solution reads p(x,y,t) = t~ 1 / 3 h(x/t 1 / 3 ), 
where h(u) = (1 — u 2 )/6 for u < 1 and vanishes for 
u > 1 Jl^,jl^Q. Notice that we recover here the i 1 / 3 
behavior observed in MD at early times for the front po- 
sition [Q. When disorder is present (go > 0), we expect 
the front to deform and eventually encounter a strong 
pinning region where it stops. The scaling of the front 
position £p with fo observed in MD simulations can be 
explained, noticing that the front will be pinned when 
the force due to the density gradient is equal to the fric- 
tion force of the order of foy/no- The gradient force can 
be estimated as aVp ~ Po/£pi where the density at the 
boundary is given by po ~ M/£ p . Thus we obtain that 

£p ~ (Ma/ foy/no) 1 / 2 = g in agreement with MD 
simulations. 

In order to confirm the validity of these considerations, 
we discretize Eq. (||) on a tilted square lattice and in- 
tegrate the equations numerically using a finite volume 
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technique with an upwind scheme to avoid numerical in- 
stabilities. At the beginning of the simulation, all flux 
lines are concentrated at the boundary of the sample (i.e. 
x = 0) and in each site of the lattice we define an uncor- 
related Gaussian pinning center density n, with mean no 
and variance iJno, thus defining a local random friction 
force. As flux lines enter into the sample, we identify 
the flux front using a burning algorithm and compute its 
average position x p . In the initial stage, x p grows as t 1 / 3 

— 1/2 

up to a crossover length scaling as £ p ~ g Q , in agree- 
ment with MD simulations (see Fig. 3). In addition, we 
measure the density profiles and find that they rescale 
with g in the same way as in MD simulations. 

The numerical integration of the diffusion equation al- 
lows for a direct analysis of the fluctuations in the front 
as a function of different internal parameters. Measur- 
ing the width W of the fronts as a function of time for 
different values of go, we find that in the initial stage W 
grows as a power law t@ where (3 ~ 0.35 until it saturates 
to a value that decreases with go. Thus the front crosses 
over from flat to fractal as it enters into the material. In 
principle we can control the strength of the fluctuations 
and the associated characteristic length £* by tuning g , 
which directly reflects experimentally measurable param- 
eters. 

In order to compare the model with experiments, we 
have to implement appropriate boundary conditions. In 
Refs. the external field was ramped at constant rate, 
which corresponds to a constant increase of the bound- 
ary density: p(0,y,t) — ht jl7|]. Integrating the equa- 
tion with this boundary condition, we observe that for 
strong disorder the flux front roughens displaying sub- 
stantial overhangs (see the inset of Fig. 4 |2J| ) . Here the 
box counting method is applied to estimate the fractal 
dimension of the front: we divide the lattice in boxes 
of size b and count the number of occupied boxes N(b), 
which decays as b~ Di , where Df is the fractal dimen- 
sion. Fig. 4 shows that for b < £* the data are well 
fitted by a power law with exponent Df = 4/3, which 
coincides with the fractal dimension of the perimeter of 
the percolation cluster. For b > £* the fractal dimension 
crosses over to Df = 1, corresponding to a flat front. The 
observed behavior is reminiscent of gradient percolation 
p6| , a model in which the concentration of occupied sites 
p decreases with the distance from the sample boundary. 
The invaded area is compact, but its perimeter is fractal 
(with Df — 4/3) up to a length scale £* which depends 
on the concentration gradient. 

In experiments on flux penetration in T^B^CuOg+z 
thin films with controlled anisotropy, Df was found to 
change continuously from Df ~ 1.5 to Df ~ 1 as the 
degree of anisotropy and the applied field was varied |j) . 
It would be interesting to check if the observed continu- 
ous variation of the fractal dimension can be explained 
by a crossover from a percolation process at short length 
scales to a flat invasion at large length scales, with a 
crossover length depending for instance on the degree of 



anisotropy. 

In conclusion, we have analyzed flux front penetration 
in disordered type II superconductors, establishing con- 
nections between microscopic flux line dynamics and a 
mesoscopic non-linear diffusion equation. This equation 
allows to interpret the results of MD simulations and pro- 
vides a useful framework to analyze experiments. 
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FIG. 1. The average position of the front, obtained from 
MD simulations injecting a constant magnetic flux from the 
boundary, is plotted as a function of time. The curve increases 
as t 1//3 and saturates at long times to a value depending on 
fo- In the inset, we show by data collapse that the pinning 
length scales as f 1,72 
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FIG. 3. The front position measured simulating the contin- 
uum equation for different values of g. In the inset we show 
the data collapse. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



f o=OM 

f o=0.02 

/„= 0.04 

/„= 0.06 

--- f o=0.10 



Pfo 





5 10 15 20 

.1/2 
X f 



150 



FIG. 2. The pinned density profiles measured in MD sim- 
ulations. In the inset we show the data collapse. 




FIG. 4. In the inset we show a density plot obtained 
simulating the continuum equation with g = 5, a — 6.28, 
no = 8, and concentration at the boundary increasing at rate 
h = 0.01. In the main figure we report the box counting plot 
of the front, averaged over twelve realizations of the disorder. 
The front has fractal dimension Df — 4/3 (solid line) and 
crosses over to Df = 1 at large length scales. 
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